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ABSTRACT 

I show in this paper that two types of magnetic torques can appear in the interaction 
between an accretion disc and a magnetic accretor. There is the well-known torque 
resulting from the difference in angular velocity between the accretion disc and the 
star, but in addition there is a torque coming from the interaction between the stellar 
magnetic field and the disc's own magnetic field. The latter form of magnetic torque 
decreases in strength more slowly with increasing radius, and will therefore dominate 
at large radii. The direction of the disc field is not determined by the difference in 
angular velocity between the star and the disc as in the Ghosh & Lamb model, but 
is rather a free parameter. The magnetic torque may therefore either spin up or spin 
down the star, and the torque changes sign if the magnetic field in the disc reverses. I 
suggest that this mechanism can explain the torque reversals that have been observed 
in some disc-fed X-ray pulsars. 
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1 INTRODUCTION 

The interaction between a magnetic star and a surrounding 
accretion disc is one of the least understood aspects of ac- 
cretion. At the same time it is of importance for understand- 
ing the spin evolution of objects as diverse as T T auri-stars 
. ^ ' and X-ray p ulsars. The generally adopted model (Ghosh & 
■ [Lamb 1979a ) has been challenged by recent observations of 
^ ' X-ray pulsars (e.g. Nelson et al. 1997). I will discuss how 
d \ the Ghosh & Lamb model fits in with current thinking re- 
garding the angular momentum transport in accretion discs 
and show that a modification of the Ghosh & Lamb model 
is necessary. 

T Tauri-stars come in two forms. Classical T Tauri- 
stars that show clear signs of having accretion discs and are 
slow rotators, and weak-lined (naked) T Tauri-stars that 
are fast rotators and lack the signature of an accretion disc. 
The weak-lined T Tauri-stars are believed to represent a 
later evolutionary stage, when the star has got rid of its 
accretion disc. Models of the rotational evolution of the T 
Tauri-stars based on the magnetic interaction between the 
disc and the star have been constructed by several groups 
(e.g Armitage & Clarke 1996, Cameron & Campbell 1993, 
Cameron, Campbell & Quaintrell 1995). 

At the other end of the range of magnetic accretors 
are the accreting neutron stars. BATSE/CGRO has made 
it possible to obtain essentially daily measurements of the 
spin frequencies of a number of the brightest X-ray pulsars. 
Some of these systems are known to possess accretion discs, 
although other objects in the group are accreting directly 



from a stellar wind. I will not pay an y atte ntion to the latter 
objects in this paper. Nelson et al. (1997) have found that 
the X-ray pulsars are oscillating between periods of spin-up 
and spin-down. The frequency derivatives are comparable 
(apart from the difference in the sign) in the two states, and 
a typical system spends about the same amount of time in 
each of the two states. The length of time a system spends 
in one of the states before switching to the other state differs 
widely between the different systems, in the case of Cen X-3 
it is 10 - 100 days, but 4U 1626-67 has only been observed 
to switch once, and thus may stay for 10 years or more in 
one state. 

The standard m odel fo r the spin changes was presented 
by Ghosh & Lamb (1979b). They describe the interaction 
of a dipolar stellar magnetic field with a diffusive disc. The 
dipolar magnetic field is penetrating the disc, which winds 
up the magnetic field in the toroidal direction, because of 
the angular velocity difference between the accretion disc 
and the star. The co-rotation radius, where the angular ve- 
locities of the disc and the star are the same, plays a key 
role in this model. The magnetic field lines penetrating the 
disc inside the co-rotation radius spin up the star, whilst 
those penetrating the accretion disc outside the co-rotation 
radius brake the star. The spin evolution of the star is there- 
fore the result of a balance between the angular momentum 
carried by the accreting matter from the disc to the star, 
the magnetic spin-up torque from the accretion disc inside 
the co-rotation radius, and the magnetic spin-down torque 
from the accretion disc outside the co-rotation radius. The 
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position of the inner edge of the accretion disc varies with 
the accretion rate such that it moves closer to the star when 
the accretion rate increases. Thus we expect the star to spin 
up, or at least spin down more slowly, when the accretion 
rate, or equivalently the luminosity, is high. 

A complicated, and not yet fully solved problem is the 
position and physical properties of the boundary layer be- 
tween the accr etion d isc and the stellar magnetosphere. 
Ghosh & Lamb (1979a) presented solutions for the boundary 
layer, but instead of solving for the magnetic field in the disc 
given a certain magnetic diffusivity they assumed a magnetic 
field a nd solved for the diffusivity. Heptinstall & Campbell 
(1998) and Brandenburg & Campbell (1995) have recently 
presented solutions for the magnetic field and velocity in the 
disc given a certain model for the magnetic diffusivity. The 
disc is terminated as the magnetic force becomes stronger 
than the viscous force in the disc (cf. Campbell 1992, 1997). 
The temperature and disc thickness increases dramatically 
close to the inner disc edge, and the disc is subject of a vis- 
cous instability similar to the Lightman-Eardley instability 
(Lightman & Eardle y 1974 Lightman 1974). 

Ghosh & Lamb (1979b) calculated the torque acting on 
the star by integrating over the surface of the disc. This is 
without a doubt the most convenient approach, but it does 
not provide much information on how the torque is trans- 
ferred to the star. To understand this problem it is necessary 
to obtain a solution for the structure of the magnetosphere. 
Simple solutions assuming no flows through the magneto- 
sphere have been presented by Bardou & Heyvaerts (1996). 
The neglect of a flow through the magnetosphere may be un - 
realistic though, and Li, Wickramasinghe & Riidiger (199C) 



have presented a class of funnel flow solutions in which the 
angular momentum is carried by the matter. Their solutions 
include a slow magnetosonic shock, in case the flow is sub- 
Alfvenic, close to the stellar surface. The shock is forced 
to be close to co-rotation with the star, which minimises 
the torque between the shock and the star. In addition the 
toroidal magnetic field at the shock is small, which keeps the 
torque between the shocked and unshocked matter small. 
The authors thus conclude that most of the angular mo- 
mentum of the matter in the fun nel h as to be propagated 
back to the accretion disc. Wang (1997) presents two major 
objections to this result. He points out that a small toroidal 
field is enough to produce a significant torque on the star. 
Furthermore if the angular momentum is propagated back 
to the disc, it must be transported outwards through the disc 
by the usual viscous torque, but he argues that the viscous 
torque is to o wea k to do this. 

Safier (1998) goes a bit further and puts into question 
whether the magnetic field of a T Tauri-star is a closed dipole 
field. He notes that the magnetic activity of a T Tauri star 
will heat up the atmosphere of the star and thus produce a 
corona and a stellar wind. This stellar wind should not be 
confused with the jets, that are typically much more massive, 
but it dominates the magnetic field of the star, which thus 
assumes an open topology as it is drawn out by the wind. 
The open field lines are mainly radial so fewer field lines 
penetrate the accretion disc than for a dipole magnetic field. 

The large uncertainties in the models make it interest- 
ing to try to simulate the interaction between the accretion 
disc and the magnetosphere numerically. Such simulations 
are in themselves sensitive to the assumptions used to con- 



struct the initial state. A critical issue is the initial magnetic 
field. Some groups have started with a stellar dipole field 
that threads the disc (Hayashi, Shibata & Matsumoto 1996, 
Goodson, Winglee & Bohm 1997), whilst other groups have 
taken into account that the disc may possess a magnetic 
field on its own (Rastatter & Neukirch 199 7, Mi ller & Stone 
1997). The simulations by Hayashi et al. (1996) and Good- 



son et al. (1997) both show that the stellar magnetic field 



is wound up by the Keplerian accretion disc and becomes 
unstable. A reconnection event follows and mass is thrown 
out more or less along the rotational axis. Goodson et al. 
suggest that this is the matter forming the stellar jet, but 
they notice that there is also a slower outflow between the 
jet and the disc, which may form a disc wind. In general the 
simulations show more evidence for outfl ow th an for accre- 
tion onto the star, though Miller & Stone ( 1997 ) find a polar 



accretion flow when their disc contains a ve rtical magnetic 
field, and likewise do Rastatter & Neukirch (1997) find that 
reconnection causes matter to be accelerated along the field 
lines leading to the polar caps. 

To summarise, the simulations are not yet mature 
enough to provide any information on the structure of steady 
accretion discs around magnetic stars. The main reason for 
this is that they have concentrated on following transient 
events for a brief time interval. However the simulations 
suggest that the existence of a dynamo-generated magnetic 
field in the accretion disc is important for the interaction be- 
tween the disc and the stellar magnetosphere. The purpose 
of this paper is to explore how the dynamo may affect the 
exchange of angular momentum between the accretion disc 
and the star. I compare the torque generated by the winding 
up of the stellar magnetic field with the torque generated by 
the coupling between a stellar magnetic field and a dynamo- 
generated magnetic field in the disc in Sect. 2. The existence 
of two different forms of magnetic torques may influence the 
spin evolution of the X-ray pulsars as suggested in Sect. 3. 
Finally my conclusions are summarised in Sect. 4. 



2 MAGNETIC TORQUES 

We can distinguish between two different ways in which the 
accretion disc and the star exchanges angular momentum; 
matter flowing over from the disc to the star carries with 
it angular momentum, and the magnetic stress at the disc 
surface exerts a torque on the disc. The torque due to the 
angular momentum carried by the accreting matter from the 
inner edge of the disc is 



M(GMr ) 1/2 . 



(1) 



where M is the accretion rate, G the gravitational constant, 
M the mass of the accreting star, ro the inner radius of the 
disc. For the rest of this paper it is useful to write 



ro 



where ta is the Alfven radius 

/ o_2..4 \ 1/7 



r A = 



2*V 



\GMM 



2 ti 



(2) 



(3) 



/io is the permeability of free space, and fi is the magnetic 
dipole moment of the accreting star. To calculate £ a de- 
tailed model of the accretion disc is needed (|Ghosh & LamtJ 
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1979a), but its value is typically around 0.5 (e.g. Frank, King 
& Raine 1992). The torque can now be written as 



Taccr — ( 27T 

where 



2\1/14 f l/2 

£ TO, 



TO = 



(GMM 2 ) fi 2 
Mo 



1/7 



(4) 



(5) 



The magnetic torque is the result of the coupling be- 
tween the vertical magnetic field of the star and the toroidal 
magnetic field in the disc. The torque acting on the upper 
surface of the disc can be written as 



2tt 



B Z B& , 

i -rdr, 

Mo 



(6) 



where B z is the vertical magnetic field, the toroidal field 
and Mo the magnetic permeability of free space. There is 
a similar contribution, but with the opposite sign from the 
lower surface, thus angular momentum is exchanged between 
the disc and the star only if B^ changes sign from the up- 
per to the lower surface. This is true if B$, is generated 
by the winding up of the stellar magnetic field, but it is 
not true for a quadrupolar magnetic field generated by a 
disc dynamo. Linear mean-field afi-dynamos with a positive 
Q-effect (e.g. Stepinski & Levy 1990, Torkelsson & Bran- 
denburg 1994a) generate preferentially quadrupolar mag- 
netic fields, but nonlinear afi-dynamos can gene rate a large 
range of different mag netic field configurations (Torkelsson 



& Branlenburg 1994b). It is not necessary that B^ is coher- 



ent across the entire accretion disc, because, as we will see 
later, the magnetic torque is concentrated to region inside 
2r . 

2.1 The torque due to the winding up of the 
stellar field in the disc 

A toroidal field is generated in the disc because the angular 
velocities of the disc and the star match only at the co- 
rotation radius, r co . Inside r co the disc rotates faster than 
the star so that B Z B$, < at the upper surface of the disc, 
and the disc is losing angular momentum to the star, whilst 
B Z B^ > outside r co and the disc is gaining angular mo- 
mentum from the star. To calculate B^ one must know the 
magnetic diffusivity in the disc (e.g. Campbell & Heptinstall 
1998, Brandenburg & Campbell 1998). A simpler approach 
is to write B^ = "/B z , where the azimuthal pitch 7 ~ 3 is 
assumed to be constant (e.g. Ghosh & Lamb 1979a). Obvi- 
ously this approach must fail close to r co , where the shear 
vanishes. Neglecting this technical complication I calculate 
the torque acting on the star due to the interaction between 
the stellar dipole field and the disc 



7M 2 2 , 
r dr - 



7M 2, 
r r dr 



Mor D 



3/iorj] 



ro 



(7) 

where fj, is the magnetic dipole moment. Substituting ro from 
Eq. (§) we get 



Tmae: — 



2 (16tt) 



1/7 



- 7 r 3 (l-2^ s 2 )ro, 



(8) 



where the fastness p arameter uj 2 = 
(Eisner & Lamb 1977) is introduced. 



{rl/GM)/{rlJGM) 



2.2 The torque due to the coupling between the 
stellar dipolar magnetic field and a 
dynamo-generated toroidal field in the disc 

It is now widely believed that the accretion is driven by a 
magnetic stress that is generated by a dynamo in the ac- 
cretion disc (e.g. Brandenburg et al. 1995, Stone et al. 1996, 
Matsumoto & Tajima 1995). The torque needed to drive the 
accretion is 



MVGMr 



(9) 



The magnetic torque, on the other hand, can be written as 

Bj,B r 



2ixr2Hr- 



where 2H is the thickness of the dis c. I write B, 
where 7d yn ~ 10 (Brandenburg 199S). Equating Eqs. 



( [i"o| ) the dynamo-generated toroidal field is 



(10) 

7dyn-Br, 

and 



B... 



/ MM07dyn } l/4 

V 4nr 2 H 



(11) 



which falls off more slowly with r than the field that is gen- 
erated by the winding up of the stellar dipole field in Sect. 



2.1 



This may be a lo wer lim it to B$ as the work by Bran- 
denburg & Campbell (1998) suggests that the torque in Eq. 
( [i"o| ) must be much larger than the torque in Eq. (|J) to re- 
distribute the angular momen tum e xchanged with the star. 

The Shakura-Sunyaev (1973) model predicts that 
H/r oc r 1//8 for a geometrically thin, optically thick disc 
in which the gas pressure dominates over the radiation pres- 
sure. Because this is sensitive to the opacity and equation 
of state, and to keep things simple, I will take H/r to be a 
constant. The torque due to the coupling between the stellar 
field and the dynamo-generated toroidal field in the disc is 
then 



Tmag.dyn 



= 2 x 2rr 



1 

r T-B<*rdr 

Mo r 1 



5 V r 



WSfHy 1 / 2 / 4TT7 dyn /i 2 



3Morg 



(GMM 2 r ) 



1/4 



of either sign. This can be re-written using Eq. O) as 



Tmag.dyn 



/r,23 4\l/28 1/2 c - 
( 2 * ) 7 d U 



V4 



-1/2 



(12) 



(13) 



This is in a sense an upper limit to the torque caused by 
the disc dynamo, as it requires the toroidal magnetic field 
to be aligned across the entire disc, but it is potentially the 
largest component of the total torque, and must therefore 
be taken into account. 



3 DISCUSSION: THE SPIN EVOLUTION OF 
DISC-ACCRETING X-RAY PULSARS 

X-ray pulsars are the best laboratories for studying the ex- 
change of angular momentum between the accretion disc 
and the accreting star, because the moment of inertia of a 
neutron star is comparatively small, and an X-ray pulsar 
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can be timed accurately. The torque is measured from the 
change in spin frequency 



27r/' 



(14) 



where / is the moment of inertia of the neutron star. Adding 
Eqs. (§), (§) and @ the spin change is 

i) = 2.4x 10~ 14 Hz s -1 I^W /7 M?/ 7 /4) 7 ?> (15) 
where 

1/2 



f = 1.2£ 1/2 +1.2^ (l-2w; 



h = J/10 40 kgm 2 , m = M/Mq, M u = M/10 14 kgs 



and ^20 = fi/10 Tm+ The Ghosh & Lamb (1979b) model 



which corresponds to the first two terms above, predicts that 
at a sufficiently low accretion rate the torque is spinning 
down the neutron star, but as the accretion rate increases 
and the inner disc edge is pushed closer to the neutron star, 
the torque changes sign, and the spin up torque increases 
with increasing accretion rate. At high accretion rates 



(17) 



It is difficult to derive any form of similar relation when the 
angular momentum exchange is dominated by the coupling 
to a dynamo-generated magnetic field in the disc. 

The amount of timing data of X-ray pulsars has in- 
creased dramati cally since the launch of BATSE/CGRO. 
Bildsten et al. (1997) have recently presented a compila- 



tion of 5 years of monitoring of X-ray pulsars. Most X-ray 
pulsars are fed by stellar winds from their supergiant com- 
panions, and are therefore irrelevant for this paper. There 
are two groups of X-ray pulsars that have accretion discs. 
These are the Be/X-ray transients and an inhomogeneous 
group of binaries with steady accretion discs. The neutron 
stars in the Be/X-ray transients sometimes pick up an accre- 
tion disc at the time of their periastron passage, which leads 
to an outburst of X-rays. The observations of transient X- 
ray pulsars, such as EXO 2030+375 (e.g. Reig & Coe 1998) 
and A0535+262 give some support to the Ghosh & Lamb 
model as the X-ray pulsars spin up during the outbursts, 
and that the spin-up rate decreases as the X-ray flux goes 
down, but a major uncertainty is the relati on between the 



obse rved X-ray flux and the accretion rate (Bildsten et al 
19970 

The persistent X-ray pulsars with accretion discs are 
comparatively few and make up an inhomogeneous set. The 
systems I will discuss in the following are 4U 1626-67, GX 
1+4, Cen X-3 and OAO 1657-415, in which the mass losing 
stars are a degenerate dwarf, a red giant, an 06-8 supergiant, 
and an OB supergiant respectively. 4U 1626-67 had been 
spinning up steadily at a rate i> — 8.5 x 1 0~ 13 Hzs" 1 from its 
discovery by Uhuru ( Giacconi et al. 1972 ) until the beginning 
of the BATSE observations, but it is now spinn ing down just 
as stea dily a t a rate v = —7.2 x 10~ 13 Hzs -1 (Chakrabarty 
et al. IS '97a ) . GX 1+4 is similar in the sense that it was also 



discovered in a state of spinning up at v — 6.0 x 10 12 Hz s 1 
(Davidsen, Malina & Bo wver 1977, Nagase et a l. 1989) , but 
since the days of Ginga ( Makishima et al. 1988 ) it has been 
spi nning d own at v = —3.7 x 10" 1 Hzs -1 . Chakrabarty et 
al. ( 1997b ) have studied the correlation between the pulsed 



X-ray flux in the 20 - 60 keV band with the spin-down rate 
of GX 1+4. There is no clear correlation, and the largest 
spin-down rates seem to occur at the highest X-ray fluxes, 
which contradicts the Ghosh & Lamb model. 

Cen X-3 shows an altogether different behaviour. On 
average it is spinning up at a rate i) = 8x 10 -13 Hzs -1 , but 
BATSE revealed a lot of fine structure in its spin evolution, 
and it is alternating between spin-up and spin-down phases 

, respec- 



tively (Bildsten et al. 1997). Typically it spends 10 to 100 



days in one state before switching to the other state faster 
than can be resolved by BATSE. OAO 1657-415 appears to 
be a sister system of Cen X-3's although the companion has 
not been identified. 

The existence of two states of comparable but oppo- 
sitely directed torques is not expected from the Ghosh & 
Lamb model. It is, though, a natural consequence, if f is 
dominated by its last term, the coupling between the stellar 
magnetic field and the dynamo-generated field in the disc. 
In that case we expect the spin-up torque to be somewhat 
larger than the spin-down torque in agreement with the ob- 
servations. 

There are two time scales that must be explained, the 
time scale over which the torque is constant, and the time 
scale for reversing the torque. These time scales are 10 - 
100 days, and less than 10 days, respectively, for Cen X-3. 
In the cases of GX 1+4 and 4U 1626-67, on the other hand, 
the torque remains constant for several years, and the sparse 
sampling at the time of torque reversal provides a generous 
upper limit for the time scale of the torque reversal. The 
time scale for reversing the torque should be the same as the 
time scale for reversing the magnetic field via diffusion in my 
model. I assume the diffusive time scale to be comparable 
to the viscous time scale 



^turb 



rr 



(18) 



where ^Wb = as sHc s is the turbulent viscosity, ass the 
Shakura-Sunyaev (1973) parameter, and c s the sound speed. 
The dynamo torque is so strongly concentrated to the inner 
part of the disc that half of the torque is due to the disc in- 
side 1.74ro. The viscous time scale at this radius is typically 
less than 0.5 days. 

It is much less clear what sets the time scale over which 
the magnetic torque is constant. The ratios of the intervals 
between torque reversals and the orbital periods may be 
roughly comparable for Cen X-3 and GX 1+4, though only 
l ower limits to the orbital pe riod of the latter are known 
( |Chakrabarty fe Roche 1997). The time interval between 



the torque reversals is comparable to the viscous time scale 
for the entire disc, which suggests that the time scale is 
set by the time over which the disc can support a given 
field configuration against diffusion. This speculation fails 
completely in the case of 4U 1626-67, which is the smallest 
of the systems I have discussed, and yet the torque remains 
constant over a time scale of several years. 

Some other mechanisms have been proposed to explain 
the transitions between spin-up an d spi n-down states in X- 
ray pulsars, van Kerkwijk et al. ( 199S ) note that the ac- 
cretion disc can be subject of a warping instability due to 
the irradiation from the neutron star. This warp may be so 
extreme that the inner part of the accretion disc flips over 
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and rotates in the opposite direction, whic h wou ld lead to 
a torque reversal. Yi, Wheeler & Vishniac ( 1997 ) suggest a 
modification of the Ghosh & Lamb model. The torque rever- 
sals are due to small changes in the accretion rate, but the 
inner part of the accretion disc changes between a standard 
thin disc and an advection-dominated flow at the same time. 
The advection-dominated flow is less efficient in radiating 
energy, which can solve the problem of the lack of a corre- 
lation between the torque and the observed X- ray fl ux. A 
third possibility was suggested by Nelson et al. ( 1997 ). The 
spin-down can be explained if the disc i s retr ograde, which 
was first suggested by Makishima et al. (1988). It is difficult 
to see how a retrograde disc can appear in a Roche-lobe 
overflow, but numerical simulations of accretion from winds 
have shown that a temporary disc rotating in the 'wrong' 
direction may appear (e.g. Fryxell & Taam 1988). 



4 SUMMARY 

I have shown that the torque acting between an accretion 
disc and an accreting star can be enhanced by the pres- 
ence of an intrinsic magnetic field in the accretion disc. The 
orientation of the magnetic field in the disc, and thus the 
direction of the magnetic torque between the disc and the 
star, is arbitrary, because the dynamo does not have any 
information on the rotation of the neutron star. In particu- 
lar it is therefore possible to reverse the torque by reversing 
the magnetic field in the disc. This mechanism may explain 
the observed torque reversals in some X-ray pulsars that are 
fed by accretion discs. The short time scale for the reversals 
compared to the long time scale over which the torque re- 
mains the same may be explained as the difference between 
the diffusive (viscous) time scales for the inner region of the 
disc, to which the torque is concentrated, and the entire ac- 
cretion disc. A problem with this connection is that it cannot 
explain the stability of the torque in the short-period X-ray 
binary 4U 1626-67. 
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